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ABSTRACT 

The angular clustering of faint field galaxies is investigated using deep imaging (/ ~ 25) obtained with 
the 10-m Keck-I telescope. The autocorrelation function is consistent with uj{0) cx d~^-^ and, although 
less steep correlation functions cannot be ruled out with high confidence, we find no compelling evidence 
for a systematic decrease in the power law index at the faintest magnitude limits. Results from a number 
of independent observational studies are combined in order to investigate the variation of the correlation 
amplitude with median /-magnitude. At /mcd ^ 23 the results obtained by different studies are all 
in rough agreement and indicate that for /mod > 22 the correlation amplitude declines far less steeply 
than would be expected from an extrapolation of the trend in the brighter samples. In particular, at 
-^mcd ^ 24 our data indicate uj{0) to be a factor ^ 7 higher than the extrapolation. A near-independence 
of magnitude is a general feature of the correlation amplitude in models in which the redshift distribution 
of the faint field population contains a substantial fraction of galaxies with z ^ 1 . In order to reproduce 
the apparent abrupt flattening of the amplitude of uj{0) observed at faint limits, approximately 50% of 
the galaxies in a sample with a depth of / 25 must be at z > 1. 

Subject headings: cosmology: large-scale structure of the universe — cosmology: observations — galaxies: 
evolution 



1. INTRODUCTION 

The angular clustering of faint field galaxies has been 
studied extensively (e.g. Efstathiou et al. 1991, hereafter 
EBKTG; Roche et al. 1993, 1996; Brainerd, Small & 
Mould 1995, hereafter BSM; Hudon & Lilly 1996; Lidman 
& Peterson 1996; Villumsen, Frcudling & da Costa 1996; 
Woods & Fahlman 1997). A prime motivation of these 
studies has been to investigate the nature of the faint field 
population. In particular, it is possible to infer the ef- 
fective correlation length of the sample and the rate at 
which clustering evolves from a combination of the ampli- 
tude of the angular autocorrelation function, uj{9), and the 
redshift distribution of the faint galaxies, N{z). These ob- 
servations can then be used to link properties of the faint 
field population with samples of local galaxies. While the 
exact interpretation remains controversial, it is generally 
accepted that overall uj{9) is fit well by a power law of 
the form d'"-^ (affhough see Infante & Pritchet (1995) 
for evidence of a fiattening in the power-law coefficient 
at faint limits). Moreover, the amplitude of uj{9) appears 
to decline strongly with apparent magnitude, although it 
remains significantly non-zero at the faintest limits. 

Some knowledge of the redshift distribution of the galax- 
ies is required in order to interpret the observed decline of 
the amplitude of lli{9) with limiting magnitude in terms 
of the evolution of real-space clustering. Although their 
precise N(z) is not known, the majority of the faint field 
population to depths of / ~ 25 appear to have (z) ~ 1 
(Kneib et al. 1996; Connolly et al. 1997). Given this con- 
straint, the rate of clustering evolution inferred from lj{9) 
is consistent with all "reasonable" theoretical possibilities. 
Similarly, the present-day correlation length of the descen- 
dents of the faint galaxy population is not constrained 
strongly. Some authors argue in favor of a value of 
close to that of local bright galaxies (e.g. Hudon & Lilly 



1995; Villumsen, Frcudling & da Costa 1996), while oth- 
ers argue for a smaller value, consistent with that seen for 
dwarf galaxies (e.g. BSM; Infante & Pritchet 1995). 

A feature which has emerged from the studies of llj{9) for 
faint galaxies in blue passbands is a slowing of the decline 
of the correlation amplitude with apparent magnitude be- 
yond B ~ 24.5-25, equivalent to a median / magnitude of 
^med ~ 22.5 (Roche et al. 1996). A change in the slope 
of the differential number counts is also seen in blue pass- 
bands at a similar magnitude (Metcalfe et al. 1995; Small 
et al. 1995). One interpretation of these features is that 
they arise from a steepening in the faint-end slope of the 
luminosity function at high redshift, z > 1 (Metcalfe et al. 
1995). Clearly it is important to confirm and extend these 
observations to test this proposal. 

Here we present an analysis oiuj{6) obtained from deep 
imaging (/ ~ 25), sufficient to probe the strength and evo- 
lution of galaxy clustering from z ~ 1 . By combining these 
observations with results from clustering analyses of other, 
brighter, /-selected samples we investigate the variation of 
uj{9) with magnitude. Using a modest extrapolation of the 
observed magnitude-redshift relation of brighter /-selected 
galaxies we then compare the observed clustering behavior 
to theoretical predictions. 

2. OBSERVED GALAXY CLUSTERING 

The data consist of deep /-band imaging of two inde- 
pendent fields (1640-^22 and 2229-^26) obtained in good 
conditions with the Low Resolution Imaging Spectrograph 
(LRIS, Oke et al. 1995) on the 10-m Keck-I telescope. 
The fields are centered on two high galactic latitude pul- 
sars and thus provide essentially random samples of faint 
field galaxies. Details of the observations, the reduction 
and cataloging of the galaxies, and the number counts of 
the galaxies are discussed in Small et al. (1995); Reid et al. 
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(1996) discuss the analysis of the faint stars in these fields. 
These frames constitute a superb dataset for the study of 
the clustering properties of faint galaxies due to both the 
excellent seeing (0.53" and 0.58" FWHM for 1640+22 and 
2229+26 respectively) and depth (50% completeness limits 
of / = 26.0 and / = 25.6 for 1640+22 and 2229+26 respec- 
tively) . To reduce stellar contamination we consider only 
objects with J > 22 in our analysis. 

The direct pair-counting method proposed by Landy & 
Szalay (1993), Loi9) = {DD - 2DR + RR)/RR, was used 
to estimate (jj{0) for all objects with 22 < I < In^. Here 
DD, DR, and RR are the number of unique data-data, 
data-random, and random-random pairs within a given 
angular separation bin centered on 0. To determine DR 
and RR, mask frames defining areas around bright stars 
and galaxies (/ < 20) were constructed for each field and 
these regions were excluded from the analysis. Also ex- 
cluded was a generous border (15 arcsec) along the frame 
boundaries, resulting in field sizes of order 30 sq. arcmin. 
Raw measurements oiLu{9) were; compiited independently 
for each field and error bars were assigned to the functions 
using bootstrap resampling of the data. Since the areas 
of the fields are almost identical and the residual stellar 
contamination is similar, it is fair to compute a mean of 
the two independent raw measurements of w(0) directly. 
The result is shown in Fig. 1 where the data points and 
error bars are the formal values obtained via a weighted 
mean of the corresponding individual functions. 

The observed number of objects is used to compute 
DD and DR. Since the area of the detector is small, 
u){6) is therefore underestimated by an amount: IC = 
f2-2/ Joj{e)dnidn2 (e.g. Croth fc Peebles 1977), the so- 
called "integral constraint" . Assuming a power law corre- 
lation function, lu{0) = Ai^0~^ with S = 0.8 as suggested 
by previous studies, we find IC1640-1-22 = (0.021 ±0.002)^^^ 
and IC2229-F26 = (0.023 ± 0.002)^^. 

Some stellar contamination remains in our catalogs and 
results in the inferred amplitude of lo{0) being lowered 
from its actual value. Fortunately the excellent seeing al- 
lowed a detailed study of the star counts in these fields 
at somewhat brighter magnitudes and from that analy- 
sis a mean stellar contamination can be estimated. Ex- 
trapolating the star counts obtained by Reid et al. (1996) 
under the assumption that beyond / ~ 22 they remain 
fairly flat (e.g. Fig. 1 of Reid et al.), we determine a stel- 
lar contamination in our catalogs of ^ 10% for 7=22.0 — 
24.0, ~9% for 7=22.0—24.5, and - 7% for 7=22.0—25.0. 
The final corrected amplitude of uj{0) is then given by 
= A]^N^y^.{Nobi - A^s)"^ where iVobj is the number of 
objects, Ns is the muiiber of stars, and A]^ is the inferred 
amplitude of uj{0) after correcting for the IC. 

The functions in Fig. 1 are consistent with power laws 
in which 6 = 0.8, as are the individual co{0) computed for 
each field. Formally the index of the best-fitting power 
law, Shest, ranges from 0.6 to 1.1 for the correlation hmc- 
tions computed from the two fields independently (Ta- 
ble 1) and from 0.7 to 0.9. for the combined fields. There 
is, however, no clear trend of Jbest with limiting magni- 
tude, nor is the formal best fit significantly better than 
that with S = 0.8. We therefore adopt oj{0) = A^0-°-^ for 
all further analysis. 
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Figure 1. Mean eorrelation funetions obtained by weighted averaging of 
the raw funetions eomputed from the two independent fields. Sohd lines 
indicate the best-fitting power law of the form lj{0) = A^0~"'^, includ- 
ing suppression due to the IC. Dotted lines show the power law that is 
formally the best fit to the observed ijj{0). See Table 1 for details of the 
individual fits. 

Table 1 lists values of uj{0) evaluated at 6 = V for the 
two individual fields and the combined sample. The er- 
rors listed are the formal errors derived from the fit of the 
power law, lo{0) oc 0~^-^, to the raw correlation functions. 
From the table, w(l') is nearly independent of limiting 
magnitude and the galaxies in 2229+26 arc slightly more 
clustered than the 1640+22 galaxies, although the ampli- 
tudes from the two fields agree at better than la. Results 
for w(l') as a function of median 7-magnitude are plot- 
ted in Fig. 2, along with those from other studies of the 
clustering of 7-selected galaxy samples. 

3. RESULTS AND INTERPRETATION 

We start by noting that the clustering amplitudes of the 
brightest samples in both our fields are in reasonable agree- 
ment with the clustering of 22 < 7 < 24 galaxies reported 
by EBKTG who find w(l') = 0.024+0.006 (corrected for 
their IC, but not stellar contamination). Allowing a mean 
stellar contamination of order 5% in the EBKTG data in- 
creases their value to w(l') = 0.027+ 0.006 (Fig. 2). 

Performing a weighted linear least squares fit to the 
data of Lidman & Peterson (1996) shown in Fig. 2, we 
find w(l') oc I^^^^ at bright magnitudes. At our faintest 
limits, however, cj(l') is ~ 7 times larger than the expec- 
tations based on an extrapolation of the linear fit to the 
Lidman & Peterson data. Moreover, for 7ined ^ 22, tj(l') 
appears to be nearly independent of magnitude, both in- 
ternally within our sample and in comparison to slightly 
brighter studies. Thus we conclude that the amplitude of 
u{0) for 7-selected samples of galaxies ceases to decline 
steeply with apparent magnitude at 7med ~ 22 and re- 
mains roughly constant at fainter magnitudes, although a 
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Table 1 
Summary of Field Statistics 



Field 


/ mag. limits 


JVobi 




u.(l')" 


X^('5bcst) '' 


x'('5 = 0.8)'' 


1640+22 


22.0—24.0 


674 


0.8 


0.038 ± 0.011 


0.22 


0.22 




22.0—24.5 


989 


0.6 


0.030 ± 0.007 


0.30 


0.39 




22.0—25.0 


1450 


0.7 


0.028 ± 0.005 


0.13 


0.18 


2229+26 


22.0—24.0 


482 


1.1 


0.050 ± 0.015 


0.03 


0.11 




22.0—24.5 


699 


0.9 


0.039 ± 0.010 


0.26 


0.31 




22.0—25.0 


1046 


0.8 


0.035 ± 0.006 


0.37 


0.39 


Combined 


22.0—24.0 


1156 


0.9 


0.043 ± 0.009 


0.13 


0.18 




22.0—24.5 


1688 


0.7 


0.034 ± 0.006 


0.31 


0.32 




22.0—25.0 


2496 


0.7 


0.031 ± 0.004 


0.28 


0.29 



a) Computed using oc ^'^ , corrected for the IC and residual stellar contamination b) P^r degree of freedom. 



modest rise or fall in the amplitude at the faintest magni- 
tudes cannot be ruled out. 

We construct a simple theoretical model to predict the 
angular clustering in our samples and compare this to the 
observations. On small scales {6-^1 rad.) an angular cor- 
relation function of the form Lij{9) = Ai^O^^ corresponds to 
a spatial correlation function of the form ^(r) = {r/ro)~^, 
where 7 = 1 + S. Observations of local bright galaxies sug- 
gest that ^(r) is fit well by a power law with index 7 ~ 1.8 
and a correlation length, ro, which varies significantly with 
both morphology and luminosity (e.g. Loveday et al. 1995 
and references therein). Following Peebles (1980), u!{9) 
and ^(r) are related through 



: r[(7-l)/2] A 
r(7/2) OT 



(1) 



where A is an amplitude factor dependent upon both 
the shape of N{z) and the evolution of ^(r) (see, e.g. 
EBKTG). The evolution of ^(r) can be parameterised by 
^{r,z) = (ro/r)'''(l-l-z)~^^"'"'^\ where e = 0.0 corresponds to 
clustering fixed in proper coordinates and e = —1.2 corre- 
sponds to clustering fixed in comoving coordinates. Recent 
measurements of ^ (r, z) for I < 22 galaxies with known 
redshifts suggest, however, that the evolution of cluster- 
ing may have been moderately rapid and that e > 0. Lc 
Fevre et al. (1996) obtain < e < 2 for the CFRS galaxies 
and Shepherd et al. (1997) obtain e 1.5 for a sample of 
CNOC field galaxies. 

A redshift distribution, N{z), is required in order to 
make a theoretical prediction of the clustering amplitude of 
the faint galaxies. As a basis for this we extrapolate from 
the N{z) observed for the CFRS galaxies (e.g. Le Fevre et 
al. 1996). For 19 < / < 22, N{z) for the CFRS galaxies 
is fit well by N{z) oc z"^ exp[— {z / zq)^] where zo is approx- 
imately equal to the median redshift. Under the assump- 
tion that the shape of N{z) does not change appreciably 
with depth, median redshifts for our faint galaxies can be 
obtained from a linear extrapolation of the CFRS Imed- 
z relation. For galaxies with 22 < 7 < /lim we expect 
z^^d = 0.81,0.86,0.91 for = 24.0,24.5,25.0 (corre- 
sponding to Jmed = 23.2,23.6,24.0). These predictions 
are similar to the limits obtained from lensing analyses of 
arclets seen through rich clusters of galaxies, which yield 
Zmcd ~ 0.8 ± 0.1 for / < 25 (Kneib et al. 1996). 

Using the above shape for N{z) with Zmed extrapolated 
from the CFRS Jmed"^ relation, theoretical predictions for 
w(l') were computed for galaxies with 19 < 7 < /um. 



Taking 7 = 1.8, Eq. (1) was evaluated at = 1' for 
e = —1.2, 0.0, 1.0 in two cosmologies: an Einstein-de Sitter 

model and an open model with fio = 0.2. The shape of the 
expected variation of log[a;(l')] with /mod is independent 
of ro (cf. Eq. (1)) and is only weakly clepcndent upon the 
cosmology. The amplitude of the function, however, is de- 
termined primarily by rg, with a small dependence upon 
the cosmology. The theoretical variation of lu(I') with Imcd 
obtained with N{z) extrapolated from the CFRS results 
in values of w(l') for a given value of /,„ed in excess of the 
linear trend suggested by bright samples, but none of the 
models reproduces the general trend of the observational 
data. As an illustrative example of this, we show in Fig. 
2 the predicted variation of a;(l') for the EdS universe us- 
ing e = —1.2 (thin solid line). In contrast to the marked 
flattening suggested by the data, beyond /med ^21 the 
predicted amplitude of w(l') decreases roughly linearly. 
Similar results are obtained for the open model and for 
e = 0.0 and 1.0. 

Due to the reduction of the proper volume element at 
high redshifts, the amplitude of uj{9) is expected to be- 
come nearly-independent of magnitude when a high pro- 
portion of the galaxies in the sample are at z > 1. Given 
that the shape of the expected variation of log[w(l')] with 
median magnitude is most strongly dependent upon N{z) 
and, in particular, upon the tail of the distribution at large 
redshift, we modify our fiducial N{z) slightly in order to 
attempt to reproduce the general trend of fiattening of the 
observed correlation amplitude at faint magnitudes. Given 
the observational uncertainties in the determination of the 
amplitude we do not attempt to constrain N{z) in a robust 
manner; rather, we simply address the question of what 
fraction of the faint field population must be at high- 2; to 
allow Lo{l') to be nearly-independent of /,ned for /med ^ 23. 
The modification consists of an extended tail of galaxies at 
z = 1-2, which shifts the median redshifts of our samples 
to z^ed = 0.92,0.96,0.98 for = 24.0,24.5,25.0. At 
-^med = 24 this places ~ 50% of the population at ^ > 1 
and ^ 17% at z > 1.5, as opposed to ~ 44% and ~ 10% 
in the extrapolation of the CFRS N{z). 

Although slight, this modification of N{z) results in a 
flattening of the amplitude of w{9) at faint limits that is 
in better agreement with the trend suggested by the data 
in Fig. 2. The theoretical curves in Fig. 2 have been nor- 
malized to reproduce the observed clustering at /mod ~ 19 
and for clarity we show only the results for the EdS uni- 
verse (but note that the results are similar for the open 
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universe). In the EdS universe, the normalizations corre- 
spond to ro = 4.0,5.0,5.9/?riMpc for e = -1.2,0.0,1.0. 
To obtain identical normalizations in the open model, the 
correlation lengths are ro = 4.4, 5.5, 6.5/i~^Mpc. Although 
none of the models provide a good fit to all of the data 
points, the general flattening trend of the data is repro- 
duced fairly well. 
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Figure 2. Observed variation of aj(l') witii median / magnitude. All 
data points have been corrected for the IC and stellar contamination. 
The thick solid line is a weighted linear least squares fit to the data of 
Lidman & Peterson (1996). Note that systematic ofTscts can be intro- 
duced into the measured values of ^^(1') by a number of effects, including 
atmospheric seeing differences between datasets, object deblending algo- 
rithms, the definition of mask regions, and stellar contamination correc- 
tions. Such offsets would typically be comparable in magnitude to the 
size of the error bars at the faintest limits. Also shown is the expected 
behavior of uj(l ) as a function of /med for a number of theoretical mod- 
els in an EdS universe. All theoretical predictions have been normalized 
to reproduce the observed clustering amplitude at /med ~ 19- The thin 
solid line shows the theoretical prediction assuming e — —1.2 and a 
faint galaxy redshift distribution extrapolated from the CFRS galaxies 
(see text). The dashed, dotted, and dot— dash lines show the theoretical 
variation of uj(l') with /med for e — —1.2, 0.0, and 1.0 using a faint 
galaxy redshift distribution that incorporates a somewhat longer high— z 
tail than the simple extrapolation of the CFRS N(z). 

4. DISCUSSION AND CONCLUSIONS 

The interpretation of the observed amplitude ofu>{9) in 
terms of a correlation length of the galaxy sample is com- 
plicated by both a lack of direct knowledge of the redshift 
distribution of these galaxies and the fact that at different 
limiting magnitudes the sample may be dominated by dif- 
ferent galaxy populations. It is entirely likely that the mor- 
phological composition of a sample will vary with depth 
due to both the shift of the restframe wavelengths sampled 
by the filter and the possibility that different galaxy popu- 
lations may evolve differently. This may contribute to the 
different conclusions drawn from clustering studies of faint 
galaxies in different passbands. At best all that can be con- 
cluded based on lo{9) analyses of deep fields is an effective 
correlation length and and an effective rate of clustering 
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evolution since it is clear that the observed clustering is in- 
extricably linked to the evolution of the galaxy population 
as a whole. Moreover, where fields are specifically chosen 
to avoid bright foreground galaxies (e.g. BSM) the inferred 
ro may be far lower than the true "universal" value since 
these lines of sight may be biased toward regions of the sky 
containing voids and, hence, would have an uncharacter- 
istically low galaxy clustering amplitude. Deep wide-field 
imaging (sufficient to average over significant amounts of 
large-scale structure) is therefore necessary to determine 
the mean amplitude of lo{0) accurately at very faint mag- 
nitudes and, hence, provide strong conclusions about the 
clustering properties of the faint field population and the 
reality of any features in the relationship of U!{9) and /med- 

Here we have measured uj{0) for faint galaxies in two 
independent random fields to / ~ 25 using high resolution 
imaging obtained with the Keck-I telescope. The angu- 
lar clustering of the galaxies is consistent ui{9) oc 9~°'^, 
though power laws which are either somewhat steeper or 
soniciwhat shallower cannot be ruled out with high con- 
fidence. Additionally, the observed clustering amplitude 
at our faint limit is consistent with that expected from 
local bright galaxies, ro ^ A-6h~^ Mpc, provided cluster- 
ing evolves in a reasonable manner and the redshift dis- 
tribution extends beyond z ^ I. This conclusion echos 
that of Woods & Fahlman (1997) from their analysis of an 
-^med ^ 23 sample. 

Fainter than /,ncd ^ 22 the amplitude of lli{9) is approx- 
imately independent of median magnitude. This result 
supports claims of a flattening in the amplitude of uj{9) 
obtained in blue samples at a comparable depth (Roche 
et al. 1996). The apparent abruptness of the flattening of 
the amplitude of lo{1') is not reproduced well by a model 
based upon a simple extrapolation of the N{z) observed 
for the CFRS galaxies. We interpret this as evidence for 
a modest increase in the fraction of the galaxies at this 
depth lying in a high-z tail compared to the extrapolation 
of the CFRS N(z). Such a population has also been used 
to explain the flattening of the differential number counts 
in the blue passbands at a depth equivalent to /med 

^22.5 

(Metcalfe et al. 1995; Small ct al. 1995). The appearance 
of this feature at the equivalent apparent magnitude in 
both blue and red passbands rules out the possibility that 
it is caused by a population of very distant galaxies, z ^ 4, 
which fall out of the sample as the Lyman-limit moves into 
the B-band. 
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